Is the direct observation of electronic coherence in electron transfer reactions possible? 
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The observabiUty of electronic coherence in electron trans- 
fer reactions is discussed. We show that under appropriate cir- 
cumstances large-amplitude oscillations can be found in the 
electronic occupation probabilities. The initial preparation of 
the system is of crucial importance for this effect, and we dis- 
cuss conditions under which experiments detecting electronic 
coherence should be feasible. The Feynman- Vernon influence 
functional formalism is extended to examine more general and 
experimentally relevant initial preparations. Analytical ex- 
pressions and path integral quantum dynamics simulations 
were developed to study the effects of various initial prepara- 
tions on the observability of electronic coherence. 



I. INTRODUCTION 

The role of quantum coherence in the dynamics of 
molecular systems in the condensed phase is still largely 
unclear. According to conventional wisdom, quantum co- 
herence should have little influence on condensed phase 
dynamics because the dephasing time of any coherent 
process is expected to be very short for systems with any 
appreciable number of degrees of freedom. The extent to 
which conventional wisdom is correct depends on many 
parameters such as the effective mass of the degree of 
freedom under consideration, its coupling strength to the 
bath, and the timescale of the bath motions. Therefore, 
the understanding of quantum coherence in condensed 
phase molecular systems can be a rather complicated is- 
sue. 

The experiments of Vos et al. on the photosynthetic 
reaction center were one of the first in which vibrational 
coherence was directly observed in an electron transfer 
(ET) system. Similar experimental observations in other 
systems have been reported since The most sur- 

prising aspect of the discovery is that quantum coher- 
ence is directly detectable and is actually preserved on 
a rather long timescale even in condensed phase systems 
with ultrafast solvent dynamics. These direct observa- 
tions of vibrational coherence in ET systems have led 
to speculations concerning the possible influence coher- 
ent quantum dynamics may have on the ET dynamics 
itself. There have been a number of theoretical studies 
on this subject P-p^. However, experiments have so 
far been unable to establish a clear connection between 
vibrational coherence and ET rates. 

In this paper, we will focus on electronic coherence 



in ET reactions. Electronic coherence, if present, would 
have a more profound effect on the ET dynamics than vi- 
brational coherence. In the absence of the environment, 
an ET system when put into a state represented by a 
superposition of two or more electronic eigenstates will 
evolve coherently. This simple picture of electronic co- 
herence becomes immensely muddled when the environ- 
ment is taken into account, because the coherent elec- 
tronic motions can be dephased by the nuclear motions 
in the environment. In fact, in most situations the elec- 
tronic motions are so strongly coupled to the nuclear mo- 
tions that the combined system is better described by the 
vibronic states. In this way the electronic motions asso- 
ciated with the electronic coherence will acquire a larger 
effective mass, and the crucial experimental question is 
whether electronic coherence can be observed directly in 
ultrafast experiments on a timescale fast enough to beat 
the dephasing time. A recent experimental attempt by 
Rcid et al. to detect electronic coherence in biruthe- 
nium mixed-valence compounds has failed to yield any 
direct evidence because the dephasing timescale seems 
to be too short even for the 20 fs time resolution of the 
experiment. Given this negative result, the possibility of 
a direct experimental detection of electronic coherence 
may appear unlikely. 

The purpose of this paper is to re-examine the theo- 
retical issues surrounding the observability of electronic 
coherence in ET systems. Contrary to prevailing skepti- 
cism, we believe that there is no intrinsic experimental 
limitation that fundamentally prevents electronic coher- 
ence from being detected. We employ simple arguments 
to show that the major difficulties encountered in the di- 
rect experimental observation of electronic coherence are 
related to how the system is prepared initially in the ex- 
periment. By using various methods of nonequilibrium 
initial preparation, the duration of time within which 
electronic coherence is preserved as well as the amplitude 
of the coherent oscillations can be increased dramatically. 
The validity of these simple ideas have been verified by 
dynamical path integral simulations. Our analyses sug- 
gest that electronic coherence can be most easily observed 
in adiabatic ET systems in which the electronic coupling 
matrix element is large compared to the frequency of typ- 
ical solvent motions. ET systems that satisfy these re- 
quirements are often found in mixed- valence compounds 
p3[ , which have typical electronic coupling matrix ele- 
ments of the order of several hundred wavenumbers. 
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The organization of the paper is as follows. Section 
H presents a simple analysis of electronic coherence in 
an ET experiment. Using this analysis, we study the 
effects of various preparation methods on the strength 
and the dephasing time of the electronic coherence in 
the adiabatic limit. Based on these results, we propose 
a number of specific ways to engineer the initial prepa- 
ration of the system so that electronic coherence could 
be more easily detected. Section [II develops these qual- 



itative ideas further by a detailed quantitative analysis 
using the two-state spin-boson Hamiltonian as a model 
for mixed-valence compounds. The Feynman- Vernon in- 
fluence functional method is extended to allow for more 
general initial preparations, and the concept of an "equiv- 
alent external field" is developed. Accurate quantum dy- 
namics simulations show that the qualitative predictions 
of the simple theories from Section |l| are generally pre- 
served when the full solvent dynamics is taken into ac- 
count. For completeness, we also include in Section IV a 
brief summary of results for the nonadiabatic region. 



II. ELECTRONIC COHERENCE FOR 
ADIABATIC ELECTRON TRANSFER 

A. Standard Preparation 

In the absence of the solvent, if the two electronic 
states before and after the charge transfer are denoted 
\L) and the electronic Hamiltonian can be written 
in the tight-binding approximation as 
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{\L){R\ + \R){L\) 
{\L)(L\^\R)(R\) , 



(1) 



where e is a bias related to the intrinsic difference in 
the redox free energies AG° of the two electronic states. 
In condensed phase ET systems, this electronic sys- 
tem is immersed in a bath. In most bimetallic mixed- 
valence compounds, the electronic coupling matrix ele- 
ments hA/2 are large enough for them to be classified as 
adiabatic ET systems . This means that the timescale 
of the pure electronic motions, set by the electronic ma- 
trix element hA/2, is much shorter than the timescale 
of the bath motions. Although the frequencies of bath 
motions often cover a broad and almost continuous spec- 
trum, there is often a typical bath frequency Wc such that 
motions with frequencies significantly higher than lOc are 
absent in the bath The adiabatic regime is then 

defined by the condition A ^ lOc- 

Under the condition A ^ cjc, we can invoke the Born- 
Oppenheimer approximation and arrive at the adiabatic 
limit of Marcus theory p^-p^. In this limit, the bath 
polarization £ can be treated as a classical coordinate, 
which couples to the dipole moment of the electronic 



state. The electronic Hamiltonian is modified by this 
coupling and it becomes 
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^+^,S]{\L){L\-\R){R\) , (2) 



where n is one-half the difference in the dipole moments 
of the two electronic states. In the adiabatic limit, the 
bath polarization is slow, so that the Born-Oppenheimer 
approximation gives as a function of the bath polariza- 
tion two electronic eigenstates with energies 
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Figure |l](a) shows the energies of the two eigenstates as 
a function of the bath polarization £ for a symmetric 
system (e = 0). 

Within Marcus theory the bath polarization obeys lin- 
ear response. Therefore, the orientation of the bath po- 
larization when the bath is uncoupled to the electronic 
system has an intrinsic gaussian distribution 
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(4) 



where /3 = l/fcsT, Zq = J d£ exp{-(3fi^£^ / X) = 
7rA//i^/3, and A is the classical reorganization energy. 
In a normal ET experiment, we imagine an initial 
preparation in which we hold the electron in state \L) and 
thermalize the bath polarization to it. In other words, 
the bath is prepared initially with a distribution 



exp[-p{fi^£yx + fi£)] 



(5) 



Related to the intrinsic bath polarization distribution (jj) , 
this distribution has the same gaussian shape, but its 
maximum is shifted to fi£ — —A/2. This is the usual ini- 
tial preparation one considers in theoretical descriptions 
of ET. Here, we shall always refer to it as the "standard 
preparation" . 

At time t = 0, we release the electron and observe its 
dynamics. Because the bath polarization is slow com- 
pared to the electron, we can compute the early-time dy- 
namics of the electron assuming that the bath is station- 
ary. This corresponds to solving for the dynamics of the 
electron using the two eigenstates for every position 
of the bath polarization £, and then averaging the result 
over the distribution function ps {£) in Eq. (H) . This pro- 
cedure is illustrated schematically in Fig. U(a). Due to 
the gaussian preparation ps, only the polarizations be- 
tween the two dashed lines are significantly populated. 
For each £, the dynamics of a biased but free two-level 
system can be solved giving PL{t]£), the transient pop- 
ulation on the L-state, with the initial condition set to 
Pl{0;£) — 1. An average over £ finally produces the 
observed Phii)'- 
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PL{t) = / d8ps{£)PL{t;£) 



(6) 



The sohd Hne in Fig. illustrates a typical result of 
Eq. for a symmetric system, e = 0. (In Fig. ||, 
time is measured in units of an inverse frequency, I/lOc, 
which for now has no special meaning other than set- 
ting a timescale.) As expected for the case of a stan- 
dard preparation, the L-state population remains ap- 
proximately unity, indicating a very slow ET rate from 
the L- to the i?-state. In Section [II A , we will verify that 
the qualitative behavior predicted by this adiabatic cal- 
culation is indeed correct using a series of exact quantum 
simulations. 



B. Nonstandard Preparations 

In this section, we will examine several experiments 
in which the initial preparation of the bath is carried 
out by a distribution function different from pg. We call 
these "nonstandard preparations" . In these experiments, 
the electron is held in the L-state for time t < as in 
the standard experiment, but the bath is prepared by a 
generalized distribution function 



Pm{£;q,f3) = ^ns exp [-/3 {p'^E^/X- q£p)] , 



(7) 



in which the bath polarization £ is equilibrated to an ef- 
fective electronic dipole moment of magnitude qp instead 
of g = — 1 as in pg. In addition, the initial bath polar- 
ization distribution can have a different width, i.e. be at 
an effective temperature /3. 

Within Marcus theory, the initial bath preparation 
specified by pNS represents the proper description for sys- 
tems which are not necessarily thermalized with the bath 
coupled to the electron held fixed in the L-state. The 
distribution in Eq. (|^) differs from that in Eq. (^) by a 
horizontal shift in fj,£ by an amount of {q + l)X/2. If 
/3 7^ /3, the distributions also have a different width. The 
following discussion focuses predominantly on the effect 
of the position of the distribution only, not its shape. 
Therefore, we will concentrate on those nonstandard dis- 
tributions with (3 — P and refer to them specifically as 
"shifted bath preparations" . In the following, we suggest 
several possible scenarios which give rise to an initial 
preparation described by Eq. (^. Similar nonstandard 
initial preparations have been considered in Refs. [ p]JT2| 
for the nonadiabatic regime, but to our knowledge non- 
standard initial preparations have not been studied for 
the adiabatic regime before. 

First, consider the experiment depicted in Fig. |l|(b). 
Using a symmetric Fe(III)-Fe(III) compound as an ex- 
ample, we imagine an experiment in which an extra elec- 
tron is photoinjected into the ion on the left at time 
t — 0. Prior to the injection, the electronic state has 
a zero dipole moment, so the bath is initially prepared 
with a distribution function pNS characterized by q = 



and (5 = (3. Using the adiabatic approximation described 
in Sect. II A , we can now compute the transient electronic 
population on the left ion PL{t) using pwo for the bath 
preparation. The result is shown in Fig. g as the short- 
dashed line. Remarkably, Ll(^) shows large-amplitude 
coherent oscillations. The origin of this coherent behav- 
ior is purely electronic, and it can be explained by refer- 
ring to the inset in Fig. 0(b). The relevant range of f , 
depicted by the region between the dashed lines, is now 
shifted to near £ = Q. Inside this range, the two adiabatic 
electronic eigenstates are close in energy (approximately 
A apart). Consequently, when forced into the L-state 
which is a superposition of the energy eigenstates, the 
electronic motion will exhibit a coherent oscillation with 
a frequency that is of the order of A. 

At this point, it is necessary to examine why preparing 
the bath thermally at q = — 1 does not seem to yield simi- 
lar coherent oscillations, whereas the nonstandard prepa- 
ration does. The absence of oscillations in the first case 
is not due to dephasing, because so far the bath has been 
assumed to be static and hence nondissipative. There 
are no oscillations, simply because for polarizations £ 
relevant for this initial preparation (provided one has 
similar parameters as in Fig. |2|), each of the electronic 
states |L) and \R) is close to one of the two energy eigen- 
states, |£^-) and |£^+), respectively, which themselves are 
of course stationary. On the other hand, if £ = 0, both 
|L) and \R) are superpositions of both energy eigenstates 
to the same extent. Therefore, when the initial prepara- 
tion provides polarizations close to zero, we can expect 
large oscillations in the time-dependent electronic occu- 
pation probabilities. 

Though elementary, this observation is the conceptual 
foundation of the present work. We see that in a usual 
experiment (i.e. with the standard preparation), lack of 
electronic oscillations is expected for reasons that have 
nothing to do with dephasing. The conventional wisdom 
that absence of electronic oscillations is due to rapid de- 
phasing is not necessarily correct, because coherent oscil- 
lations can be suppressed by the way the system is pre- 
pared instead of by dephasing. Moreover, the standard 
way ET experiments are carried out are not the proper 
experiments for observing electronic coherence. We want 
to re-examine this issue in detail. From our analysis of 
the adiabatic limit we can understand the physical pic- 
ture and try to construct experiments which are best po- 
sitioned to detect electronic coherence directly. In Sec- 
tion [II we go beyond the adiabatic limit and formulate 



the problem rigorously including the quantum dynam- 
ics of the bath as well as shifts in the initial prepara- 
tions. With that formulation, we will verify that even for 
moderately large A/wc, the adiabatic predictions remain 
qualitatively correct. 

We can imagine a few more experimental situations in 
which nonstandard initial preparations arise. In the ex- 
periments by Vos et al. [Q , for example, the donor state is 
populated by a fast laser pulse which leaves the environ- 
ment essentially unchanged during excitation. Therefore, 
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the initial bath polarization distribution corresponds to 
a gaussian wavepacket that is equilibrated with respect 
to an electronic state different from \L). Such an initial 
preparation can be characterized by a pNS with s.q ^ —1. 
In Fig. 1^, we have also shown PL(t) for a nonstandard 
preparation with q = —0.5. This could be achieved, say, 
by ultrafast laser pulse excitation from a ground state 
with an electronic dipole moment of magnitude — 0.5/x. 
Even for this situation, we find that Phit) (long-dashed 
curve in Fig. ||) still displays some, albeit much weaker, 
coherent oscillations. 

Another experiment in which a nonstandard bath 
preparation can be achieved is an ultrafast back electron 
transfer (b-ET) experiment depicted in Fig. ^c) . In this 
experiment, a fast laser pulse is used to excite a mixed- 
valence compound into the charge-transfer (CT) band to 
initiate a fast b-ET reaction back to the ground state. 
When the system is placed on the CT state, the polar- 
ization will quickly move away from the excitation region 
and may not have sufficient time to equilibrate. This sit- 
uation is describable by a nonstandard bath preparation 
Pns with q = +1. 

The few experiments described above may not always 
be practically feasible. For example, in the experiment 
illustrated in Fig 0(a) , only in a highly idealized situation 
will the photoinjected electron accquire a pure \L) state 
instead of a superposition of \L) and \R). There is no 
doubt that the proper initial preparation will be tricky 
to establish, otherwise electronic coherence would have 
been seen already. Alternative ways to achieve similar 
shifted pre parat ions will need to be formulated. 

In Sect. [II B , it is shown that remarkably the same 
ET behavior associated with shifted preparations of the 
type of pNS in Eq- @j can also be achieved by using the 
standard preparation, but coupling the electron to a suit- 
ably chosen external time-dependent electric field that is 
tailored to mimic the effects of the shifted bath. For ex- 
ample, instead of releasing the electron into a polarized 
environment, one could imagine an initially unpolarized 
bath and apply an external field that couples to the elec- 
tric dipole moment just as the bath polarization would, 
thus producing identical ET behavior. As the effects of 
the initial bath preparation dissipates, we would have 
to adjust the external field to mimic the time-dependent 
effects of the bath. We call such a field the "equiva- 
lent external field" , because it exactly compensates for 
the effects on the ET that a different initial preparation 
would have had. This method is depicted schematically 
in Fig. |l|(d). A time-dependent electric field like the one 
depicted in Fig. |l](d) is difficult to generate experimen- 
tally, because the decay time of the field must be similar 
to the bath relaxation time which is often of the order of 
a few hundred femtoseconds. An electric field of sufficient 
strength that varies on that timescale can probably only 
be generated optically. We will show in Section III, that 
approximating the true equivalent electric field crudely 
by a static field plus a CW laser field and applying it 
to the equilibrium system can induce electronic coherent 



behaviors very similar to that with a g = initial prepa- 
ration. Such a field is schematically shown in Fig. |l|(e). 

Finally, there may be additional difficulties with the 
experimental observability of electronic coherence. Even 
if the proper initial preparations or equivalent electric 
fields can be achieved, the coherent oscillations in the 
electronic populations may be difficult to detect experi- 
mentally. Since the coherent oscialltaions, if present, are 
expected to have a frequency ~ A, given the best time 
resolution of currently available lasers of approximately 
20 fs, this limits the systems in which coherent ET dy- 
namics can be observed to those having A smaller than 
a few hundred wavenumbers. 



III. THEORY: PATH INTEGRAL 
FORMULATION AND QUANTUM MONTE 
CARLO SIMULATIONS 

A. Nonstandard Initial Preparations 

In this section, we want to verify that the qualita- 
tive predictions from the last section concerning the de- 
tectability of electronic coherence are still correct when 
the full quantum dynamics of the bath are taken into 
account. 

ET with a quantum mechanical bath has thoroughly 
been studied theoretically in the framework of the spin- 
boson model | p^ , p^j2l| . In the spin-boson model, the 
electronic part of the Hamiltonian is described by Eq. (||), 
but the bath polarization /if = CkXk is repre- 
sented by a linear combination of an infinite collection 
of linearly- responding effective solvent modes {xk} [ p^ . 
An infinite number of solvent modes is required for a 
proper description of a truly dissipative environment. 
The coupling constants Cfc that determine the interaction 
strength of the solvent modes with the electronic states 
are specified by a spectral density J{uj) PH . 

Within this model, the ET dynamics, that is the 
time dependent population probabilities Pl(0 {PR{t)) 
of the left (right) state can be determined numerically 
exactly using a method known as the quantum Monte 
Carlo (QMC) path integral method. This technique has 
been successfully employed to investigate the dynamics 
in many ET systems and has been described at length 
in many of our previous papers, e.g. Refs. [^3|-^5|, and 
reviews p6| , ^ . We will not repeat the details here. 

We imagine an ET experiment in which the electron 
is created in the donor state prior to the start of the 
reaction. If the electron is held in the donor state for 
a period of time long enough for the bath to equilibrate 
to it, the initial density matrix is given by a 'factorized" 
form (l|j2l|: 



p{0) = \L){L\ 



(8) 



where Hb is the bare bath Hamiltonian. In typical situa- 
tions, the electron is strongly solvated in the donor state. 
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and the initial preparation specified by Eq. represents 
the standard way an ET experiment is often carried out 
pGf . In the classical limit, this is equiv alent to the stan- 
dard preparation we considered in Sect. [I A and depicted 
in Fig. |l|(a). 

Nonstandard initial conditions result when the bath 
does not have sufhcient time to equilibrate to the electron 
created in the donor state. An example of an experiment 
in which this may occ ur is the one depicted in Fig. |l](b) 
and considered in Sect. [IB . Following the photoinjection 
of an electron into the donor state in Fig. l|(b), the bath 
may not have enough time to equilibrate to the electron 
if the subsequent ET rate is rapid. This would result in a 
factorized initial condition similar to Eq. (^), except the 
bath is equilibrated to a zero-dipole-momcnt electronic 
system such that the preparation is now defined by the 
initial density matrix 



p(0) = \L){L\ e 



-PHe 



(9) 



This gives rise to the nonstanda rd sh ifted bath prepara- 
tion we have considered in Sect. [IB. In general, we can 
represent any factorized system preparation by an initial 
density matrix in the form 



piO;q) = \L){L\ e 



-P{HB-qti£) 



(10) 



where the parameter q determines whether the bath 
preparation is standard or shifted. Obviously, g = — 1 
corresponds to the standard preparation and g = to 
the experiment of Fig. |l|(b). The same p(0) can also be 
used to describe the initial preparation in the b-ET ex- 
periment in Fig. |l|(c), the appropriate value of q there 
being +1. For a laser pulse excitation from the ground- 
state Born-Oppenheimer surface, at least in principle, 
any value of q can be generated. 

The computer programs from the previous QMC stud- 
ies 1 33 -3^] were modified to include the nonstandard 
preparation defined by Eq. (|l^) . A series of QMC simula- 
tions were performed for a number of ET parameters and 
different nonstandard preparations to examine the influ- 
ence of shifted bath preparations on the time-dependent 
occupation probability Ph^t)- In particular, we looked 
for conditions under which large-amplitude oscillations 
were evident in Phit)- 

Since the bath obeys linear response, all effects of the 
bath on the electronic states can be described by the bath 
polarization correlation function, which is related to the 
spectral density J{lij) by 

,c-^. Nc-^, ^^ ^ , r, ,COS'h(Lu\h(3/2~it]) 

^ ^ ^ ^ " (2^)27r Jo smh{ujh[3/2) 

(11) 

For most ET systems, a sufficiently realistic choice for 
J(w) is the ohmic spectral density with an exponential 
cutoff, 



J{uj) = {2Trha)Lue 



j/u. 



(12) 



where the dimensionless friction parameter a is related 
to the classical reorganization energy A by A = 2ahuJc. A 
solvent with an ohmic spectral density has a bath polar- 
ization correlation function that is approximately con- 
stant at short times and (at low temperatures) decays 
algebraically at long times. This behavior is qualita- 
tively similar to what has been observed in experiments 
for several solvents and in molecular dynamics sim- 
ulations [p3|-p5[. A bath with an ohmic spectral density 
describes a polarization that by itself is overdamped and 
consequently does not intrinsically exhibit vibrational co- 
herence. This ensures that oscillations in the electronic 
populations whenever detected must be due solely to elec- 
tronic coherence, but not due to vibrational coherence of 
the nuclear coordinates. 

It is noteworthy that the same ohmic-like spectral den- 
sity has been found by simulations to characterize the 
primary ET in bacterial photosynthesis Therefore, 
regarding the primary charge separation ET process, vi- 
brational coherence is expected also to have no relevance. 

We shall now compare QMC data to the predictions 
from Sect. || to verify that the general conclusions 
reached there remain correct when a fully dynamical and 
dissipative bath is included in the model. First, we ex- 
amine QMC results for a symmetric (e = 0) ET sys- 
tem with standard and shifted bath preparations. All 
frequency and energy parameters can be expressed in 
dimensionless units of lOc, and for this calculation, we 
have selected A/lOc = 4, (3hujc = 3 and a friction con- 
stant a such that the classical reorganization energy is 
X/TiLOc — 20. The parameters are representative of many 
mixed- valence compounds, which generally have lJc of the 
order of a hundred to several hundred cm^^. The QMC 
results are shown in Fig. ^. For a typical lJc — 100 cm"-'^, 
one LOct on the time-axis corresponds to 53 fs in real time. 

The solid line in Fig. ||(a) shows QMC data for the 
standard bath preparation q = —1. As expected, be- 
cause the temperature and the reorganization energy cho- 
sen put the system in the activated region, the ET rate 
is slow and the time-dependent occupation probability 
Phii) decays so slowly that on the timescalc plotted it 
appears to remain almost constant except for a fast in itial 
transient. The qualitative prediction made in Sect. HA 
for the same set of parameters is also shown in the same 
figure as the dashed line. Obviously, the adiabatic the- 
ory prediction captures most of the qualitative features 
of the early-time dynamics of the ET. 

Figure ^(c) shows QMC results for a nonstandard bath 
preparation with q — Q. This corresponds to a bath 
initially equilibrated to a zero-dipole-moment electronic 
state, such as the experiment depicted in Fig. 1(b). In 
this case, there are large-amplitude oscillations in the 
electronic population, indicating clearly the presence of 
electronic coherence. The close resemblance of the ex- 
act QMC results to the qualitative adiabatic theory pre- 
diction for the same parameters (dashed line) indicates 
that the reasoning given in Sect. p^I B| for the detectability 
of electronic coherence in experiments with nonstandard 
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bath preparations is indeed correct. Figure |^(b) shows 
results for an intermediate value of g = —0.5. Coherent 
oscillations are still present, though to a smaller extent. 

The discrepancies between QMC data and the adia- 
batic theory predictions visible in Fig. 3 are as expected. 
Since the adiabatic limit does not have dissipation, the 
electronic coherence persists indefinitely there, whereas 
the oscillations observed in the QMC data disappear 
rather quickly. Consequently, an exact QMC calculation 
is necessary to estimate whether the coherence predicted 
by qualitative adiabatic theory is indeed preserved for a 
long enough period of time so that oscillations can be 
detected in experiments. Another discrepancy between 
QMC and adiabatic prediction is that the oscillation fre- 
quency in the QMC results is initially similar to the adi- 
abatic prediction, but then increases with time. This 
behavior is also easily rationalized with the picture of a 
dynamical polarization field which moves away from the 
barrier region with time into regions where the differ- 
ence between the energy eigenvalues is larger, resulting 
in higher oscillation frequencies. 

We have shown results for only one set of parameters. 
A large number of QMC calculations on other param- 
eter sets have also been performed, and we found that 
similar coherence effects can be observed for many other 
systems as well, as long as A is similar to or larger than 
bJc- Our goal is to highlight the possible electronic coher- 
ence effects that could be detected in experiments, but 
not to give a detailed range of parameters within which 
electronic coherence can be observed (such a task would 
have been overwhelming). Therefore, we will not explic- 
itly show results from the other calculations here. 

Results presented so far apply to symmetric systems. 
The qualitative behavior can change quite dramatically 
depending on the bias e. Figure ^ shows effects of stan- 
dard and nonstandard bath preparation in the activa- 
tionless region where fit = A. In the activationless case, 
if the bath was completely static, even a standard prepa- 
ration would place the bath inside a region where the 
two energy eigenstates are close to each other, and adi- 
abatic theory would therefore predict large coherent os- 
cillations. Clearly the QMC data in Fig. ^ for g = — 1 
show no such oscillations. Thus, adiabatic theory delivers 
a qualitatively incorrect picture for this situation. The 
reason why adiabatic theory is wrong is not difficult to 
understand. When the bath is placed right at the acti- 
vationless crossing point on the Marcus parabola where 
there is no thermal barrier to cross from the donor to 
the acceptor state, the bath should move rapidly to form 
the product. The adiabatic assumption that the bath is 
slow in the activationless region (and regions close to it) 
is incorrect, rendering the qualitative prediction of adia- 
batic theory invalid. In reality, the bath being dynamical 
moves quickly away from the crossing region into regions 
where the difference between the energy eigenvalues is 
much larger, thus destroying any sign of coherence. 

Next, we consider QMC results for the b-ET experi- 
ment depicted in Fig. |l|(c). The physical situation is the 



following. For all t < 0, the electron is held fixed in the 
\R) state. Then, at i = 0, the electron is excited to the 
\L) state by an ultrafast laser pulse, such that the posi- 
tion of the environment remains essentially unchanged, 
i.e. with the nonstandard preparation characterized by 
q — +\. Figure |^ shows results for a b-ET experiment 
with the same parameters as in Fig. ^, except a nonstan- 
dard bath preparation oi q — -t-1 was used. Once again, 
coherent electronic oscillations are clearly observable. Al- 
though the electronic coherence is evident in Fig. ||, its 
origin is more subtle than the results exhibited in Fig. ^ 
and this deserves a closer examination. 

When the system is excited to the \L) state, the bath 
being dynamical will move down the Marcus parabola to- 
ward the crossing region. For this reason, the predictions 
of adiabatic theory, which assumes the bath is completely 
static, are invalid just as in the activationless case. To un- 
derstand the origin of the coherence, a dynamical picture 
is required. We focus on the sequence of events following 
the vertical excitation of the system into the \L) state. 
As the bath moves toward the crossing region, the elec- 
tron being coupled to the bath will also evolve. There 
are two different scenarios. If the bath moves slowly, the 
electron will have time to adjust to the bath's position 
and quickly attain a mixture of the CT state and the 
ground electronic state to approach one of the two en- 
ergy eigenstates. In this case, the system would exhibit 
no electronic coherence, because the electron is already 
in an energy eigenstate. On the other hand, if the bath 
moves rapidly, the electron may not have time to adjust 
and would stay in the \L) state, which is not an energy 
eigenstate. After a brief time delay, which is evident in 
the QMC data in Fig. ||, the bath polarization arrives at 
the crossing region where the two electronic energy eigen- 
state are closest in energy, and the electronic population 
would then exhibit large-amplitude coherent oscillations. 
In other words, whether electronic coherent oscillations 
are detectable in the electronic occupation probability 
depends crucially on the electronic coupling A and the 
bath dynamics timescale l/wc- We are in the best posi- 
tion to observe electronic coherence in b-ET experiments 
if A is larger than or comparable to lOc- 

B. The Equivalent Electric Field 

Within the spin-boson model, the path integral formu- 
lation easily reveals that the effects of the shifted initial 
preparations on the ET equal those of an appropriately 
chosen time-dependent external field. The reasoning goes 
as follows; When considering an initial preparation as in 
EQ. ([lo|), the time evolution of the density matrix to a 
time t is just a special case of the expressions given in 
in that the electronic variable is kept fixed at 
value q throughout the imaginary time path. This leads 
to an influence functional, in which the shift of the prepa- 
ration only appears in a term i{l + q) (It R{t)x(t) . 
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Here, X = Q ~ q' q, q' are the electronic variables 
on the forward and backward paths respectively. R is 
the imaginary part of a function Q{z), which is (2ii/h)'^ 
times double integral of the correlation function Eq. (|ll| ) 
with Q(0) = 0. 
In particular. 



(13) 



cosh(n/3w/2) - cosh(a;[(n/3/2) - iz] 
sinh(S,/3u;/2) 



Integrating by parts and noting that x(0) = x(0 = [the 
electron is assumed to be in the \L) state initially, and the 
population in state \L) is measured at time t], this term 
becomes + q) dTR{T)x{T). For a strict ohmic 

spectral density (ujc — > oo), which is a good approxi- 
mation for dissipative environments most often encoun- 
tered in solid-state systems, the integral vanishes and not 
surprisingly the preparation has no effect. This is why 
this kind of term is usually tacitly omitted in connection 
with NIBA calculations ||l^j2^. We will see, however, 
that for parameters that are typical for many chemical 
systems, this term does have a significant effect on the 
ET, even in regimes where the NIBA is still valid, see 
Sect. IV. From the Marcus picture it is clear that the 
prefactor, 1 + q, is proportional to the displacement of 
the polarization fr om t he standard preparation that we 
defined in Section II A. If q = — 1, the bath is equili- 



brated to the donor state \L), i.e. the bath is prepared 
in the standard fashion. In this situation the term above 
vanishes exactly. It is very advantageous for computa- 
tions that the influence functional, although reflecting 
certain initial preparations, does not include integrations 
over imaginary-time paths. 

For a physical understanding, it is helpful to note 
that a preparation with q gives each path a weight 
exp[i(l -I- q) /q dTR{T)xiT)]. On the other hand, the ac- 
tion of the bare TLS includes a factor exp[z J* dr ex{T)]. 
Obviously, (1 -I- q)R{T) acts just like a time-dependent 
bias that could be due to a strong electric field coupHng 
to the dipole moment associated with the electronic co- 
ordinate. Studying the ET with a nonstandard initial 
preparation is therefore manifestly equivalent to study- 
ing the ET of the same system with a standard initial 
bath preparation (q = —1), but with a time- dependent 
external field that mimics the preparation effects. We 
call this field the "equivalent electric field" . The physical 
reason for this observation is that for the ET the source 
of the field coupling to the electronic dipole moment is 
unimportant, be it due to the polarization of the solvent 
or to the external field or both. 

The specific form of the equivalent external field can 
be seen by noting that the phase associated with each 
electronic path is unchanged under the transformation 



(e,g)^(e + (l + g-)i?(T),-l) 



(14) 



Therefore, the equivalent external field appears as an ad- 
ditional term to the Hamiltonian, 

Ai? = |(1 + q)R{t) (\L){L\ \R){R\^ . (15) 

Eq. (|l^) establishes the connection between nonstandard 
initial preparations and the work already done on driven 
spin-boson systems, see Refs. p^-|4l|. It is interesting to 
note that Eq. ( |l5| ) decreases in time just like the classical 
damping kernel 7(t) since from Eq. ( p^ ) 

7(<) = (2/7r) / dw^^cos(wi) = 2hRit) . (16) 
Jo ^ 

So far, the equivalent external field has just been a 
name for the effect that nonstandard distributions have 
on the influence functional. However, as mentioned in 
Sect. [IB, the equivalent external field may actually be 
applied to the system to mimic a nonstandard initial 
preparation in an experiment in which the actual prepa- 
ration is the standard one, as in Figs. 0(d) and (e). To 
give a flavor, we want to repeat the calculation in Fig. ^ 
for the standard preparation but under the influence of 
an equivalent external field. The proper form of the field 
is given by Eqs. (13) and ( p^ ) and shown in the inset 
of Fig. ^ as the dashed line. In a real experiment, it is 
much easier to generate a sinusoidal field. Therefore, in 
the QMC calculation, we used a simple cosine field plus 
a static bias shown as the solid line in the inset of Fig. ||. 
With this approximate equivalent electric field, we see 
from the solid line in Fig. ^ that the external field induces 
a behavior very similar to the q = initial preparation. 

Instead of applying the cosine external field one might 
try something even simpler and couple a CW source into 
the system where the phase at t = is selected at ran- 
dom. The long-dashed curve in Fig. ^depicts the average 
of the electronic population from four different calcula- 
tions, in which the external fields all had a cosine time- 
dependent modulation but had phases shifted by 0, 1/4, 
1/2 and 3/4 of a period. Although much weaker, coher- 
ent oscillations are still clearly evident. A typical field 
strength required would be about W^KV/cm, which is 
not difficult to achieve. 

What we have shown up to this point is that given spe- 
cific initial preparations, electronic coherence can lead to 
oscillations in the electronic populations that are pro- 
nounced enough to be detected in experiments. How- 
ever, a femtosecond time-resolved spectroscopic analysis 
of the electronic populations is not a trivial experiment 
and the interpretation of the signal is often complicated 
by effects unrelated to the ET itself. Moreover, in some 
systems where we would expect the ET to be coherent, 
the oscillations are on the timescale of tens of fs, and 
therefore cannot be resolved even in the fastest pump- 
probe experiments possible today. Finally, we want to 
point out that our calculations are based on the assump- 
tion of a gaussian bath, neglecting e.g. the role of strong 
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anharmonic modes that couple to the ET. Clearly, na- 
ture is more complicated than our model, but it is safe 
to conclude that strong dissipation by itself does not pre- 
vent electronic coherence from being observed. In spite 
of the above limitations, we hope that our results will en- 
courrage experimental attempts to detect electronic co- 
herence in the scenarios we have suggested. 



IV. INITIAL PREPARATION EFFECTS FOR 
NONADIABATIC ELECTRON TRANSFER 

In this section, we briefly discuss the effect of non- 
standard initial preparations in the nonadiabatic limit, 
A ujc- The effects should be smaller in this limit, be- 
cause as argued above when the timescale for a electron 
hop, 1/A, is much longer than the typical timescale of 
of the bath, l/ujc, the initial bath state is not expected 
to matter. Studying the effects of shifted initial prepa- 
rations in the nonadiabatic limit becomes especially in- 
teresting in the context of vibrational coherence, because 
a polarization distribution shifted away from the energy 
surface minimum may lead to an interesting dynamical 
behavior. 

In the nonadiabatic limit, the ET transfer can be 
very accurately described in the framework of the non- 
interacting blip approximation (NIB A) [Q. Using the 
concept of the equivalent electric field, a nonstandard 
preparation can be easily incorporated in the NIBA. In 
fact, the proper form of the NIBA expression for PL{t) 
in the presence of a driving field has already been given, 
e.g., in Refs. |3|,^, though from a very different per- 
spective. All we need to do is to modify these expressions 
by choosing the interaction energy AH{t) according to 
the external field in the form (15), getting 



PLit) = -(AV2) / dre 



-Sir) 



sin(i?(r))x 



sm er 



-f (1 + q)[R{t) - R{t - t)] + cos(i?(T)) X 



cos er + (1 + q)[R{t) - R{t - r)] (2Pl(< - t) - 1 



(17) 



We first examine the short time behavior of the ET. 
For times smaller than l/coc, Ri^) is proportional to 
the reorganization energy, A, and at high temperatures, 
S{t) = Xr'^/n^p. With PL{t) « 1 for t < l/w^, one finds 



^^g-ArV?i^/3cos([(A/?i)g-f e]r) 



(18) 

The "rate" T{t) defined here displays a few oscillations 
of frequency [{\/h)q + e] at extremely short times. This 
frequency corresponds to the vertical distance of the di- 
abatic Marcus parabolas at = gA/2, so these oscil- 
lations are due to electronic resonance. In the nonadia- 
batic limit, the ET takes place on a much longer timescale 



compared to this resonance. As a result, although these 
oscillations can be observed in Phit), they exist only for 
very short times (up to w 1/wc) and their amplitudes are 
only 1% of Phit)- Therefore, electronic coherence plays 
only a minor role in nonadiabatic ET. 

The more interesting question is whether a nonstan- 
dard initial bath distribution causes effects on timescales 
t « l/wc. Experiments in the photosynthetic reaction 
center have given evidence of vibrational coherence on 
the same timescale as the ET itself , where oscillations 
of the nuclear coordinate distribution (a "wavepacket" ) 
within the donor energy parabola could be observed be- 
fore dephasing occurred. Since the nuclear coordinate 
£ must lie in the Landau- Zener crossing region for ET 
to occur, the nuclear motions should result in oscillatory 
transfer rates and in the electronic occupation probabil- 
ities. As already pointed out, with an ohmic bath there 
is no vibrational coherence. However, one should still 
be able to detect the relaxation of the polarization from 
the transfer rate. To test this wavepacket picture, we 
chose to investigate a rather small electronic coupling, 
A — QAujc, with A = "ihtOc and h(3wc = 10. 

Figure shows Phit) for g = 0, — 1 and —2 obtained 
from a NIBA calculation (solid lines), and from quantum 
dynamics simulations (dashed lines), q = — \ corresponds 
to standard prepartion with the bath equilibrated to the 
donor state, g = corresponds to a polarization distribu- 
tion initially placed at the crossing point of the Marcus 
parabolas (inset of Fig. |^) where the majority of the ET 
occurs. The ET proceeds rapidly. Finally, q — —2 cor- 
responds to a polarization distribution centered on the 
opposite side of the minimum. To get to the crossing re- 
gion, the bath has to first relax to the parabola minimum, 
causing a delay in the ET. This is most easily discerned 
from the rate r(<). Figure ^ shows the difference of the 
rates for q — —2 and from the standard preparation 
q — —1. They are clearly out of phase with each other, 
in agreement with the picture of two wavepackets start- 
ing at opposite turning points. 

A different situation is found in the activationless 
regime, where the Landau-Zener crossing region is now 
at the bottom of the donor parabola. Two wavepackets 
starting at opposite turning points should result in rates 
that are in phase. Figure |9| shows that this is indeed the 
case for times up to l/wc- Once the polarization reaches 
the crossing region, the two situations are no longer sym- 
metric. The inset of Fig. || shows a closeup of the cross- 
ing region. Compared to the one approaching from the 
right (5 = 4), the polarization coming from the left side 
(g — —6) will have an easier time tunneling, by simply 
sliding down the lower energy surface. This is reflected 
by the transfer rates shown in Fig. ^. We see that in 
the nonadiabatic regime nonstandard preparations also 
have a visible effect on the ET consistent with a sim- 
ple wavepacket picture. However, the effect is much less 
significant than in the adiabatic regime. 

Whether vibrational coherence causes oscillations in 
the populations or not depends on the specific spectral 
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density that is involved. If the spectral density is domi- 
nated by a narrow range of frequencies, not surprisingly 
we find oscillations in PLit) when a nonstandard initial 
preparation is used. These oscillations vanish when the 
polarization distribution is initially placed at the mini- 
mum of the donor parabola. The effect of vibrational 
coherence on the ET can be predicted by looking at the 
relaxation behavior of the initially shifted polarization 
distribution, i.e. the form of R{t) or the damping ker- 
nel. For the photosynthetic reaction center, the classical 
correlation functions have been determined by molecu- 
lar dynamics simulations . Although they show some 
high frequency oscillatory modulations, the basic behav- 
ior is an overdamped relaxation. We conclude that vibra- 
tional coherence should not cause significant oscillations 
in the electron populations in the reaction center, even 
if the initial wavepacket is largely displaced. This is also 
the conclusion from quantum dynamics simulations em- 
ploying the spectral density determined in |24|], and in 
agreement with recent experimental findings while 
another group suggests the opposite ji^ . 

Not surprisingly, taking a model spectral density con- 
taining a large concentration of weights within a very 
narrow range of frequencies, our NIBA and QMC calcu- 
lations both showed that PL{t) becomes oscillatory. 

V. CONCLUDING REMARKS 

The unambiguous detectability of electronic coherence 
in ET reactions is still an open question from an exper- 
imental point of view. In this paper, we have shown 
that there are situations, in which one might be able 
to indeed observe oscillatory behaviors in electronic oc- 
cupation probabilities which are caused by electronic 
coherence. The most appropriate systems seem to be 
symmetric mixed-valence compounds, which are close 
to the adiabatic ET regime. The crucial point is then 
to consider preparations where the solvent polarization 
initially forms a distribution centered near the Landau- 
Zener crossing region. We have given some examples for 
how one can create such an initial preparation in practice. 
Before the wavepacket slides away, i.e. before dephasing 
occurs, the electron can have enough time to oscillate 
back and forth several times between the electronic sur- 
faces. The resulting large-amplitude oscillations should 
be easily detectable in experiments. 

Besides the experimental relevance, we believe that 
these results might be of use for other theoretical stud- 
ies of nonequilibrium preparations in condensed phase 
systems. We have extended the Feynman- Vernon influ- 
ence functional method to account for nonstandard initial 
preparations. 

From the path- integral expressions, one can easily de- 
duce that such nonstandard preparations can be thought 
of as an equivalent external field. This establishes a link 
with the field of driven dissipative quantum systems and 



this connection is useful in actual computations. An- 
other important consequence of this concept is that one 
can quite easily decide whether a particular ET system 
may exhibit oscillatory rates by inspecting the equiva- 
lent external field. If the latter does not show oscilla- 
tions but simply relaxes to zero, the polarization motion 
will be overdamped, and oscillatory behaviors must have 
an electronic origin. For that reason, we conclude that 
vibrational coherence cannot cause oscillations in the oc- 
cupation probabilities for a large class of common ET 
systems. A prominent example is the primary ET step 
in the photosynthetic reaction center. 

To conclude, we hope that this paper will stimulate 
experiments investigating electronic coherence in electron 
transfer reactions, as well as theoretical studies dealing 
with the effects of special initial preparations, e.g. due to 
the laser pulse excitation of a wavepacket, in condensed 
phase systems. 
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FIG. 1. (a) The electronic energy eigenvalues E±{£) as 
a function of the bath polarization £. The shaded gaussian 
represents ps{£) that arises from a parabolic free energy sur- 
face (Marcus parabola). The mixed- valence compound for 
this standard preparation is sketched. The excess electron 
sits on the left metal center, and the bath polarization has 
adjusted to it (indicated by the density of dots), (b) Exper- 
imental situation where we expect large electronic coherence 
effects. For all times t < 0, there is no electric dipole moment 
and therefore {£) = 0, depicted by the uniform distribution of 
dots. At time t = an electron is injected into the left atom 
putting the electron in state \L), but with pi<is{£;q = 0,f3) 
(see inset), (c) The initial preparation for backward ET in 
which the bath is equilibrated with respect to the excess elec- 
tron in the acceptor, (d) While the actual system is prepared 
in the standard way, the initial bath polarization can be com- 
pensated by a suitable equivalent external field, so that the 
ET behaves just as in the case of Fig. |l|(b) . (e) In the actual 
experiment, it will be necessary to emulate the bath prepa- 
ration using an optically produced equivalent external field 
atop a static background. 



FIG. 2. Phit) evaluated in the adiabatic limit, for stan- 
dard bath preparation ps{£) (solid curve) and shifted bath 
preparations pNs(f;g ~ —0.5,/? = /?) (long-dashed curve) 
and pisis{£',q — 0, P — P) (short-dashed curve). The param- 
eters used were A = 4cJc, A = 20hu)c and hPuc — 3 for a 
symmetric system, e = 0. 
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FIG. 4. QMC results for Phit) in the activationless regime 
he = A. Other parameters are the same as in Fig. El 



FIG. 3. Comparison of the data from Fig. ^ (dashed lines) 
with the corresponding QMC results (solid lines). 
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FIG. 5. QMC data for Pi(t) in the case of a b-ET reac- 
tion. Parameters are the same as in Fig. H. 
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FIG. 6. Pl for the standard initial preparation but with 
an equivalent external field. In order to reproduce ET of 
the g = initial preparation (short-dashed curve given here 
for comparison), the external field must take the form of the 
dashed line in the inset. If the equivalent field is approxi- 
mated by a sinusoidal field on a static background (solid line 
in inset), the ET (solid line) behaves very similarly to the 
nonstandard preparation showing large electronic coherence. 
The long-dashed line is an estimate of the population if a 
random-phase cosine-modulated external field is used instead 
(see text). 



FIG. 7. PL{t) for A = 0.4wc, A = Shuc and hPuic = 10 
for initial bath preparations with q = 0,-1,-2 in the sym- 
metric case € = 0. The solid curves have been obtained in a 
NIBA calculation, the short-dashed curves by QMC. The in- 
set shows the position of the corresponding initial polarization 
distributions in the Marcus picture. 
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FIG. 8. The normalized transfer rates r'^{t) for the situa- 
tion of Fig. 0. 




FIG. 9. The NIBA-calculated normalized transfer rates 
in the activationless regime fit = X. The wavepackets start 
at opposite sides of the donator parabola minimum. The in- 
set shows a closeup look of the Landau- Zener crossing re- 
gion, which allows to explain the asymmetric behavior in the 
wavepacket picture. 
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